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Abstract 



We present a unified treatment of the quantum mechanics of 5-factory 
and neutrino oscillation experiments. While our approach obtains the 
usual phenomenological predictions for these experiments, it does so with- 
out having to invoke perplexing Einstein-Podolsky-Rosen correlations or 
non-intuitive kinematical assumptions. 
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1 Introduction 



The quantum mechanics of mixing during propagation is at the heart of both B- 
factory and neutrino oscillation experiments. In this paper, we will treat both these 
experiments in the same way. Our treatment has several advantages. In dealing 
with the S-factory experiments, it avoids having to invoke real but nonetheless puz- 
zling Einstein-Podolsky-Rosen correlations. In dealing with neutrino oscillation, our 
approach avoids the non-intuitive assumption that all the interfering neutrino mass 
eigenstates in a beam have the same energy. 

2 B-Factory Experiments 

In a typical 5-factory experiment studying CP violation, an electron and positron 
collide and form an T(4s). The T(4s), a bb bound state with intrinsic spin S = 1, 
then decays to a pair of B mesons. We are interested in the case where these B 
mesons are a B° and a B°, neutral particles that then undergo B° — B° mixing. 
Viewing the process in the T(4s) rest frame, as in Fig. (TJ we are especially interested 
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Figure 1: One of two coherent contributions to the process in which an T(4s) decays 
into a neutral B pair, after which one B decays semileptonically to a negative lepton 
l~ and other particles X, and the other B decays to a hadronic CP eigenstate fcp- 
The various features and factors appearing in the figure are defined in the text. 



in events in which one of the B mesons decays semileptonically into a negatively 
charged lepton t~ plus other particles X at a spacetime point (ti,xi) in the T(4s) 
rest frame, while the other B meson decays hadronically into a CP eigenstate fcp at 
some other spacetime point (t C p,xcp) in the T(4s) rest frame. We shall calculate 
the amplitude for the entire sequence, from the T(4s) decay through the B decays 
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of interest pp. In doing so, we shall make use of the fact that if an unstable particle 
has mass m and width T, adding up to a complex mass A = m — iT/2, then the 
amplitude for this particle to propagate for a proper time r in its own rest frame is 
exp(— iXr) pp. (If this propagation is through a distance x during a time t in some 
frame in which the particle has momentum p and energy E, then the non-decaying 
part of exp(— 2 At), namely exp(— imr), is simply the familiar quantum-mechanical 
plane- wave factor exp[i(px — Et)].) 

In order to employ the simple propagation amplitude exp(— iXr), we work in the 
B mass eigenstate basis. There are two neutral B mass eigenstates: the heavier one 
Bh, and the lighter one Bl- Since these B mesons are spinless, while the T(4s) has 
5=1, the decay T(4s) — > BB leaves the B mesons in a p wave. Since one cannot 
have two identical spinless bosons in a p wave, one of the daughter B mesons must 
be a Bh, while the other is a Bl- There are then two scenarios that contribute to 
the sequence T(4s) — > BB — > (£~X)( y f C p). In the first, pictured in Fig. [JJ it is a B H 
that heads for the point (te,xi) and there decays into £~X, and a B L that heads for 
the point (tcp,xcp) and there decays into fcp- In the second scenario, it is a Bl 
that travels to (te,xe) and decays into £~X, and a Bh that travels to {tcp,xcp) and 
decays into fcp- Since the mass difference mp H — mp L = 3 x 10~ 4 eV is exceedingly 
tiny, one will not know which scenario was responsible for any given event. Therefore, 
the amplitudes for the two scenarios must be added coherently. 

For the scenario in which Bh yields £~X and Bl yields fcp, the amplitude is the 
product of factors shown in Fig. [TJ The first of these factors, A(B H to £; B L to CP), 
is the amplitude for T(4s) to decay into a Bh headed towards the point (t£,Xi) and 
a B L headed towards (tcp,xcp)- The only feature of this amplitude that is relevant 
for us is that, owing to the antisymmetry of the p wave, it is antisymmetric under 
Bh <=>■ B L - Next, there is the amplitude exp(— iXhTi) for the B H , with complex 
mass Xh, to travel from the T(4s) decay point, which we shall call (0,0), to the 
point (t£,xe). The quantity ti is the proper time that elapses in the Bh rest frame 
during this journey. There is also the analogous amplitude exp(— iXlTcp) for the B L , 
with complex mass A^, to travel from (0,0) to (t C p,xcp)- The quantity r C p is the 
proper time that elapses in the Bl rest frame during this trip. Finally, there are the 
amplitudes A(B H — > £~X) and A{B L fcp) f° r the B H to decay to £~X and the 
Bl to fcp- Thus, including the scenario not pictured in Fig. [fl the amplitude Amp 
for T(4s) — > BB, following which one B — > £~X after traveling for a proper time t#, 
and the other B — >■ fcp after traveling for a proper time tqp, is given by 



Amp oce~ iXHTe e~ iXLTCP A(B H -> e~X)A(B L -> f C p) - (B H B L ) . (1) 



We note that this Amp is Lorentz invariant. 

Using the fact that Bh and Bl have essentially the same width T, we may write 
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where m = (rng^ +mg L ) / 2, and Am e (tyi b h — Tnp L ) . For the Bh,l decay amplitudes, 
we use the standard relation 
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(3) 



and 



in which Sq^m * s a CP-violating B° — B° mixing phase coming from the CKM quark 
mixing matrix. Since only a B°, but not a B°, can decay semileptonically to a 
negatively- charged lepton £~, this relation immediately yields the relevant factors in 
A(B HiL — > £~X). For the decays to fcp, we assume that B° — > f C p is dominated by 
a single diagram, so that we may write the B° and B° decay amplitudes as 

A(B° f C p) = Me i5 cKM e iasT (4) 

A(W f C p) = r] f Me- iS cKM e ia sT _ ( 5 ) 

Here, M is the magnitude of the dominating diagram, Sq KM is its CP-violating CKM 
phase, asT is its CP-even strong interaction phase, and i]f is the CP parity of the CP 
eigenstate f C p. 

From Eqs. CQ) - ©, we find that the probability P[T ^ BB ^ (£-X){f CP )} that 
T(4s) — > BB, and then one B — > t~X after and the other B — > fcp after r^p, is 
given by 

P[T ^BB ^ (rX)(f CP )} oc |^lmj9| 2 oc e - r(rcp+T ^{l sin 0sin[Am(r CP - r^)]} . 

(6) 

Here, <j) is the CP-violating phase defined by 

<P^^ckm + ^ckm) ■ (7) 

The probability of Eq. fl6]) is the usual result j2], except that times in the T(4s) rest 
frame are replaced by proper times in the B rest frames. This is a negligible correc- 
tion. Much more importantly, we have derived the usual result without invoking any 
puzzling Einstein- Podolsky- Rosen correlations. In the usual treatment, one imagines 
to begin with that the decay B — > £~X, at a time ti in the T(4s) rest frame, is the 
first of the two B decays to occur. Then, since only a B°, but not a B°, can undergo 
the decay B — > £~X, and since one cannot have two identical spinless bosons in a 
p wave, one concludes that at the time ti in the T(4s) rest frame, the surviving B 
must be a pure B°. This surviving B then undergoes B° — B° mixing, and evolves 
according to the Schrodinger time evolution equation, until it decays into fcp- This 
is a perfectly legitimate approach, but it does leave one wondering how the surviving 
B "knows" how and when the other B, with which it is not communicating, decayed. 
Our approach avoids this puzzle. 
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3 Neutrino Oscillation 



Now let us treat neutrino oscillation in the same way as we have treated 5-factory 
experiments. Consider the illustrative process pictured in Fig. |2j There, a neutrino 
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Figure 2: The contribution of neutrino mass eigenstate vi to the process in which a 
pion undergoes the decay tt — > fi + u, after which the muon interacts with matter 
at some spacetime point, and the neutrino interacts in a detector at some other 
spacetime point, producing an electron e and other particles X. The various features 
and factors appearing in the figure are defined in the text. 



is produced via the two-body decay tt — > // + v. We shall describe the process in the 
pion rest frame. The muon travels from the spacetime point where the pion decays, 
which we shall call (0,0), to the spacetime point (t^x^), where it interacts with the 
matter surrounding the pion decay region. The neutrino journeys from (0, 0) to the 
spacetime point (t u ,x u ), where it interacts in a detector and creates an electron. The 
propagating neutrino is one or another of the neutrino mass eigenstates z/j. 

As we did for the S-factory experiments, we shall calculate the amplitude for the 
entire sequence shown in Fig. [2j For the decay tt — fi + z/j, the relevant factor in the 
amplitude is U* { , where U is the Pontecorvo-Maki-Nakagawa-Sakata leptonic mixing 
matrix — the leptonic analogue of the CKM quark mixing matrix. Similarly, for the 
creation of an electron in the detector by z/j, the relevant factor in the amplitude 
is U e i. For the interaction of the muon in matter, we shall call the amplitude S^. 
For the propagation of the neutrino z^, the amplitude is exp(— im^r^), as explained 
when we discussed B meson propagation. Here, m % v is the mass of h>i, and r* is 
the proper time that elapses in the z/j rest frame during its journey from (0, 0) to 
(t u ,x u ). In the phase m l u T l v1 we are neglecting the extremely tiny neutrino decay 
width. However, we are allowing, through the index i on r*, for the fact that the 
proper time that elapses in the vi rest frame while the Ui travels from (0, 0) to the given 
interaction point (t u , x v ) depends on the z/j energy in the pion rest frame, hence on 
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the Vi mass, and consequently on which Vi is involved. [In the S-factory experiments, 
the analogous effect — the dependence of T£ or r C p on whether it is a B H or a B L 
that is propagating — is negligible. For example, in the Bh propagation amplitude 
exp(— iXnTt), the factor Xh = m + Am/2 — iT/2 certainly involves a term of 1st order 
in Am = rriB H — rriB L , but the dependence of the proper time Ti on Am is only 2nd 
order in (m_e H ~ m B L )/{ m B H + m B L ) — 3 x 10~ 14 .] Finally, for the propagation of the 
muon, the amplitude is exp(— iX^r^), where A M = — iY^/2 includes both the muon 
mass m M and its width T^, and is the proper time that elapses in the muon rest 
frame during its journey from (0, 0) to the point (t^, x M ) where it interacts. By giving 

an index i, we are taking into account the fact that the muon and the neutrino 
are kinematically entangled. The energies of both of these particles in the pion rest 
frame depend on the mass of the emitted neutrino, so that they depend on which 
it is. Consequently, the proper times that elapse in the rest frames of the muon and 
the neutrino both depend on which is emitted, at least in principle. 

Since we cannot make measurements that determine which Vi was involved in 
any given event without destroying the oscillation pattern, the amplitudes for the 
contributions of the different Vi must be added coherently. Thus, including all the 
factors shown in Fig. [2J the amplitude Amp for the entire sequence of events pictured 
there is given by 

Amp= £ ^^K-^r^^-imX^ (g) 

i=l,2,3 

Here, we assume for the sake of illustration that there are just three neutrino mass 
eigenstates z^. We note that the amplitude of Eq. ([8]) is Lorentz invariant. 

How do the muon and neutrino propagation amplitudes exp[— i{m^ — iT fl /2)r^\ 
and exp(— im^r*) actually depend on the index i! In the propagation amplitude for 
the muon, 

r;=— (E^-^x,) , (9) 

where and are, respectively, the muon energy and momentum in the pion rest 
frame when the emitted neutrino is z/j. In evaluating this expression, we choose a 
muon interaction spacetime point x^) at which x M and are related by 

^ = v% = ^t, . (10) 

Here, and E® are, respectively, the velocity, momentum, and energy the muon 

would have in the pion rest frame if all neutrinos were massless. The spacetime points 
passed by the peak of the wave packet describing the muon propagation in greater 
detail than is needed here would satisfy Eq. f fTUj) to an excellent approximation. 
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From Eq. fl9]), the difference between the muon travel proper times corresponding 
to 7r — > n + Vi and n — > jj, + Uj is given by 



T — T 



(11) 



Now, for given neutrino mass m„, the energy E^ of the recoiling muon in the pion 
rest frame is given by 



E„ 



m 2 + m 2 — m\ 



M 2m % 
where m v is the pion mass. Thus, we have 
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Thus, to lowest order in the squares of the neutrino masses, 
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From Eqs. (TTTj) . ( flOj) . (|T3|) . and (fl5|) . we then have, to lowest (i.e., first) order in the 
squares of the neutrino masses, 
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(16) 



That is, to lowest order, the muon propagation amplitude, exp(— im^rp exp[— (r^/2)r*], 
actually does not depend on i. This was first noticed by Akhmedov and Smirnov [3]. 
The factor exp(— zm^r*) will have no effect on the absolute square of the amplitude of 
Eq. (|8]) for the sequence of events in Fig. |2j The factor exp[— (r At /2)r*] ) which decays 
with time, will influence the dependence of the event rate for the sequence on how 
long the muon travels before it interacts, but it will not affect the neutrino oscillation 
pattern. Since it is only the latter in which we are ultimately interested, we can drop 
the entire muon propagation amplitude from the amplitude of Eq. OH]). 

In the propagation amplitude for the neutrino v^, exp(— zm^r*), we have 



E v t u p v x v 



(17) 



where E % v and p l u are, respectively, the i/j energy and momentum in the pion rest frame. 
Since in practice neutrinos are ultra-relativistic, we choose a neutrino interaction 
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spacetime point (t u , x v ) at which t v = x v = L°. From the relation E % u = [m 2 + 
[m l u ) 2 — m 2 l ]/2m 7T , we find that the energies of two different neutrino mass eigenstates 
Pi and Uj differ by 

^- B ' = i Am « • (18 > 

From the relation (p l u ) 2 = {E l u ) 2 — (ml) 2 , we find that, to lowest order in Am|, the 
momenta of V{ and Uj differ by 

Here, E® is the energy the neutrino would have in the pion rest frame if it were 
massless. From Eq. ( FlTl) - ( TT9l) . it then follows that the propagation phases of two 
different mass eigenstates and Uj differ by 

mi4-mi4 = [(Ei-El)-(pi-pl)}L 

L° 



Ami 



i + — ^ 

E?, 



2m 7T 

Thus, we may take the neutrino propagation amplitude exp(— im^rl) to be 



e l{m ^ 2 ^ v ; (21) 



and all the relative phases in our amplitude of Eq. (jSJ) for the process of Fig. [2] will 
be correct. 

Using this result, and dropping the muon interaction and propagation amplitudes, 
which, being i- independent, will not affect the neutrino oscillation pattern, we find 
from Eq. (ED that 

Amp= Yl U;,e HmJ WU* . (22) 

1=1,2,3 

We recall that the L° and E® appearing in this expression are the neutrino travel 
distance and energy (neglecting neutrino mass) in the pion rest frame. However, if 
neutrino oscillation vs. travel distance in the laboratory frame is to be observed, then 
obviously the location in the lab of the neutrino source — the pion in the process of 
Fig. |2| — must be known to within an oscillation length or so. But then, by the uncer- 
tainty principle ApAx > %, there must be some spread in lab-frame pion momenta 
[I]. The pions whose decays yield the neutrinos of an oscillation experiment cannot 
all be at rest. However, because our neutrinos are ultra-relativistic, the travel time 
of one of them in the pion rest frame is essentially equal to its travel distance L°. 
Similarly, its momentum in the pion rest frame is equal to its energy E®. Then, from 



7 



the Lorentz transformation, the lab-frame travel distance L and energy E of this 
neutrino are related to their pion-rest-frame counterparts, L° and E®, by 



L ^ ov(L° + /^L ) _ L° 
E 7w (£° + p v E°) El 



(23) 



Here, (3 n is the velocity of the parent pion in the lab, and 77r = 1/yl — PI- From 
Eq. (123]) . we see that in terms of lab-frame variables, the Lorentz-invariant amplitude 
of Eq. (J22D is 

Amp= U^e-^yAuei . (24) 



i=l,2,3 



The sequence of events depicted in Fig. [2] is what would usually be called — > v e 
oscillation, with the addition of an interaction between matter and the muon that 
recoils against the neutrino in the pion decay. We are interested primarily in the 
probability for the — > v e oscillation, integrated over all the possible fates of the 
muon. Apart from a possible overall normalization factor, this muon-integrated — > 
v e oscillation probability, Piy^ — > z/ e ), will be given by the absolute square of the 
amplitude of Eq. ( 124"|) . from which the muon interaction and propagation amplitudes 
have been removed. Assuming that the mixing matrix U is unitary, we find from 
Eq. (J23J that 



P{v^v e ) = -4^K(f/;^ e ^^)sin 2 (AmJ- 



4 ^ 

+2 J2 9 (U; t U et U M U: 3 ) sin(Am2. A) . (25 ) 

This result agrees completely, even in normalization, with the usual expression for 
P{ v n ~^ v e) In this paper, we have derived -P(^ M — > u e ) in the same way that we 
treated 5-factory experiments. We allowed for the kinematical entanglement between 
the muon and the neutrino produced in a pion decay [6], but this entanglement proved 
to be irrelevant. In deriving Piv^ — > u e ), we did not need to make any assumption 
about how the energies of the different coherent mass eigenstate components of a 
neutrino beam are related. Instead, the energies of our coherently contributing mass 
eigenstates followed from energy-momentum conservation in tt — > \w. 



4 Summary 

In summary, we have treated 5-factory and neutrino-oscillation experiments, two 
very interesting examples of the consequences of coherence in quantum mechanics, 
within a common framework. In the processes studied in 5-factory experiments, 
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puzzling Einstein-Podolsky- Rosen correlations really do exist, but in our approach 
one can understand these processes without directly encountering those perplexing 
correlations. In treating neutrino oscillation, our approach allows for quantum en- 
tanglement but yields the standard result for the oscillation probability. It arrives 
at this result without having to make any non-intuitive assumption about neutrino 
mass-eigenstate energies. 

The approach described in this paper raises some interesting issues that will be 
explored in future work. 
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